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Abstract.  In  this  paper  we  study  a  time-harmonic  fluid-solid  interaction  model  problem 
in  one  dimension.  This  is  a  Helmholtz-type  system  equipped  with  boundary  and  transmis¬ 
sion  conditions.  We  show  the  existence  of  a  unique  solution  to  this  problem  and  study  its 
stability  and  regularity  properties.  We  analyze  the  convergence  of  finite  element  methods 
with  respect  to  appropriate  energy  norms.  Computational  results  are  also  presented. 


1.  Introduction 

We  shall  analyze  a  model  problem  for  a  fluid  solid  interaction  proccess  in  one  di¬ 
mension.  Here,  we  assume  that  we  are  given  a  layered  fluid-solid-fluid  medium  with 
configuration  =  [0,  L],  L  >  0.  12  is  divided  in  three  subintervals  ft  =  fii  U  O2  U  123 
where  fti  =  ^^3  =  [^2^L]  will  be  the  ‘‘fluid”  parts  and  122  =  [^ij^2]  will 

be  the  “solid”  part.  We  assume  that  the  time-harmonic  acoustic  wave  equations  are 
satisfied  in  the  fluid  and  the  time-harmonic  elastic  wave  equations  in  the  solid: 

Pxx  +  k^P  =  -9i,  in  fti 

(1.1)  (aui)x  -f-  Pgu  =  -/,  in  0,2 

Pxx  +  k^P=-92,  infis- 


Here  k,  k  R,k  >  ko  >  0  denotes  the  wave  number  and  p  =  p{x),p  =  p{x)  denote  the 
pessure  in  Oi  and  Oz  respectively  and  u  =  u(x)  the  displacement  in  O2.  We  assume  that 
a  and  g  are  positive  constants  that  characterize  the  materials,  and  gi,g2if  3.re  given 
functions.  We  assume  radiation  conditions  at  the  boundary  of  O 

Pxi0)  +  ikp{0)  =  0, 

ML)  -  ikm  =  0. 


Also  the  continuity  requirement  for  the  pressure  and  for  the  displacement  yields  the 
transmission  conditions  at  Xi  : 


(1.2b) 


Pi(xi)  -  k'^u{xi)  -  0, 

p(xi)  +  aux(xi)  =  0, 
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and  at  X2  : 


^1(^2)  -  k^u{x2)  =  0, 
p{x2)  -\- aUx{x2)  =  ^  ■ 

(1.1-2)  is  a  system  of  three  equations  in  one  dimension  equipped  with  boundary  and 
transmission  conditions.  Each  one  of  these  is  a  Helmholtz-type  equation  where  the 
wave  number  k  appears  as  a  common  parameter  of  the  system.  We  shall  propose  a 
weak  formulation,  show  the  existence  of  a  unique  solution  for  this  problem  and  study 
its  stabihty  and  regularity  properties.  Further  we  consider  finite  element  methods  for  the 
approximation  of  the  solution  and  analyze  the  convergence  behavior  of  these  methods 
with  respect  to  appropriate  energy  norms.  Throughout  this  work  a  particular  emphassis 
is  given  in  the  dependence  of  the  constants  on  the  frequency  parameter  k. 

We  shall  state  a  weak  formulation  of  this  problem  below.  First,  we  introduce  some 
notation:  Let  £f’"(5),  m  positive  integer,  be  the  usual  Sobolev  space  of  complex  valued 
functions  which  have  derivatives  in  the  sence  of  distributions  up  to  the  order  m.  Let 
II  •  II  and  (•,  •)  be  the  L^(5)-norm  and  inner  product,  respectively.  We  shall  denote  by 
H  the  space  x  H^{n2)  x  If  q,v  and  q  are  three  smooth  test  functions 

defined  on  ^1,^2  and  fls,  we  see  that  (1.1)  and  (1.2)  give 

(1.3)  /  Pxq^dx  -k^  I  pqdx  -  k^u{xi)q{xi)  -  ikp(0)q{0)  =  /  giqdx  , 


(1.4)  /  aUxVx 

Jcii 


dx  —  k^  I  guvdx  —  p{xi)v{xi)  +  p{x2)v{x2)  =  / 
JU2  dUi 


fvdx , 


(1.5)  /  pxqxdx-k^  pqdx  +  Pu{x2)q{x2)-ikp{L)q{L)  =  g2qdx . 

jQa  dOa 

Multiplying  (1.4)  by  k"^  and  summing  (1.3),  (1.4)  and  (1.5)  we  arrive  to  the  following 
weak  formulation  of  problem  (1.1-2):  Find  lA  =  (p,  u,^  G  'hi  such  that 

(1.6)  H(W,V)  =  (:F,V)o,w,  VV6W 

where  for  V  =  {q,v,^  £  'H,  ^ 

B{U,V)=  I  pxqxdx-k^  I  pqdx-k^u{xi)q{xi)-ikp{0)q(0) 

Jfii  Jiii 

(16a)  J  guvdx  -  k^p{xi)v{xi)  +  Pp{x2)v{x2) 

+  /  Pxlxdx-k^  /  pqdx  +  k^u{x2)q{x2)  -  ikp{L)q{L) 


00 


and 


(1.6b)  (^,V)o,w  =  (51,9) +  + (5^2,1)  • 

In  the  sequel  we  shall  see  that  (1.6)  is  indeed  a  suitable  variational  formulation  for 
the  problem  at  hand.  We  show  that  problem  (1.6)  admits  a  unique  solution,  we  study 
its  stability  properties  and  we  analyze  the  convergence  of  numerical  approximations  to 
it.  In  order  to  describe  the  results  of  this  work  we  shall  need  some  more  notation.  For 
k,  k  >  0,  fixed  smdU  =  {p,u,p),  V  =  (q,  u,?)  €  let 

(1.7)  {U,V)o,'H  =  iP,Q)  +  k^{u,v)  +  {p,q). 

Of  course  (•,  •)o,n  defines  an  inner  product  in  :=  £2(^2)  x  L^{ilz)  equivalent 

to  the  standard  inner  product  of  this  space.  If  V  is  an  element  of  'H ,  by  Vi  we  denote 
{Qxj'^xtQx)  where  q,  v  and  q  are  the  components  of  V  .  Further  by  (•,  we  shall  denote 
the  inner  product  in  7i  defined  by 

(1.8)  (Z7,  V)i,w  =  {Ux,Vx)o,'H  +  ^^(^5  V)o,-H  . 

Finally  ||  •  ||o,-h  and  ||  •  lli.w  shah  be  the  norms  induced  by  the  inner  products  (•,  ■)o,n  and 
(-,  Oi.W  •  These  norms  are  better  suited  to  the  nature  of  our  problem  and  furthermore 
help  us  to  reduce  the  dependence  on  k  of  the  constants  in  the  stability  estimates. 

In  Section  2  we  study  the  continuous  problem  (1.6).  Using  variational  techniques, 
we  show  that  if  .F  €  L^,  and  U  is  a.  solution  of  (1.6),  then 


(1.9)  ||Z^||i,W  <  C'i|l.F||o,w  , 

where  Ci  is  a  constant  independent  of  k.  The  estimate  (1.9)  establishes  the  uniqueness 
of  the  solution  of  (1.6)  and,  combined  with  the  fact  that  B(-,  •)  satisfies  a  Garding-type 
inequality,  yields  the  existence  of  2A  as  weh,  by  a  standard  alternative  argument  of  the 
theory  of  differential  equations.  Using  an  estimate  of  the  type  (1.9)  for  a  properly  chosen 
auxiliary  problem  we  prove  the  Babuska  -  Brezzi  (BB)  condition  for  the  bilinear  form 

B{-)  : 


(1.10) 


KtB{U,V) 

V6W 


1 

-^k 


1,-w 


where  7  is  a  positive  constant  independent  of  k.  Furthermore,  using  (1.9)  and  (1.10) 
we  derive  other  similar  “regularity”  estimates  that  are  very  useful  in  the  convergence 
analysis  of  the  numerical  methods  for  our  problem. 

In  Section  3  we  consider  Galerkin-finite  element  methods  for  approximating  the  so¬ 
lution  of  (1.1-2),  If  Sh  is  a  suitable  finite  element  subspace  of  H  consisting  of  piecewise 
polynomial  functions,  we  define  the  Galerkin  approximation  ZY/i  G  Sh  to  lA  as  the  solution 
of 


(1.11) 


=  W£Sh. 
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We  establish  existence  of  a  unique  solution  of  (1.11)  provided  that  the  quantity  —k 
is  small  enough,  where  h  is  the  maximum  of  the  mesh  sizes  and  sq  is  the  polynomial 
degree  of  the  functions  of  Sh  ?  cf.  Section  3.  Further,  under  the  same  hypothesis  on  /i, 
we  show  that  the  discrete  analog  of  the  BB  condition  (1.10)  is  satisfied  on  Sh  a^nd  a 
quasioptimal  estimate  of  the  form 

(1.12)  \\U ml 

is  valid,  where  C»  is  a  positive  constant  independent  of  h  and  h.  Thus,  since  the  approx¬ 
imation  properties  of  Sh  will  be  known,  (1.12)  will  imply  convergence  of  optimal  order 
of  the  finite  element  approximations  in  the  ||  •  ||i,'H  norm.  This  result  is  independent  of 
the  choice  of  the  particular  method  and  of  the  finite  element  spaces  {h  version,  p  version 
OT  h  —  p  version).  Note  that  stability  estimates  similar  to  (1.9),  derived  in  section  2  are 
needed  in  the  analysis  of  the  discrete  problem  and,  in  particular,  in  the  proof  of  (1.12). 

In  Section  4  we  present  numerical  results  on  the  behavior  of  the  finite  element  approx¬ 
imations  to  (1.1-2)  (coupled  problem).  As  it  is  expected  from  the  theoretical  analysis, 
increasing  the  polynomial  degree  sq,  one  decreases  the  influence  of  the  pollution  effect 
on  the  error  and,  of  course,  increases  the  rate  of  convergence.  These  results  are  com¬ 
pared  to  analogous  computational  results  obtained  previously  for  the  single  reduced 
wave  equation  (uncoupled  problem),  [IB  1,  2],  [BS].  We  also  investigate  whether  the 
restrictive  condition  is  small”,  mentioned  above,  is  actually  needed  or  not.  The 

experiments  of  Section  4  provide  numerical  evidence  of  the  fact  (for  So  =  1)  that  the 
restriction  hk^  <  a,  a  constant,  is  necessary  for  the  quasioptimal  estimate  (1.12)  to 
hold.  On  the  other  hand,  the  results  of  [IB1,2]  for  the  single  reduced  wave  equation 
indicate  that  this  restriction  is  rather  pessimistic  for  establishing  existence-uniqueness 
and  stability  (BB  condition)  of  the  finite  element  approximations.  Indeed,  for  the  un¬ 
coupled  problem  and  in  the  case  where  Sh  consists  of  piecewise  polynomial  functions  in 
a  uniform  grid,  it  is  shown  in  [IB  1,2]  that  the  discrete  BB  constant  is  proportional  to 
provided  that  kh  remains  bounded. 

Work  on  Galerkin  approximations  for  a  model  problem  of  fluid-solid  interaction  has 
been  also  done  by  Demkowicz,  [Dj.  A  general  variational  setting  was  considered,  for 
which  the  asymptotic  convergence  of  the  Galerkin  approximations  was  estabhshed.  Also, 
in  a  numerical  evaluation  of  a  one  dimensional  coupled  problem  it  was  observed  that  the 
discrete  (and  continuous)  inf-sup  constant  decreases  algebraically  with  the  wave  number 
/j.  Particular  attention  is  given  in  [D]  to  the  sensitivity  of  the  stability  constant  with 
respect  to  wave  damping  at  the  interface  of  the  two  media.  The  approach  we  take  here 
is  different  and  focusses  mainly  on  the  stability  of  the  continuous  problem  (1.1-2)  and 
on  the  stability  and  convergence  of  the  discrete  problem  (1.11)  with  respect  to  the  wave 
number  k.  One  may  say  that,  for  the  coupled  problem  and  without  special  assumptions 
on  the  finite  element  method,  our  results  are  an  extension  of  similar  results  known  for 
the  one  dimensional  single  reduced  wave  equation,  cf.  eg.  [BGT],  [AKS],  [DS],  [IB  1,2]. 
Note  however  that  we  do  not  use  explicit  representations  ofU  and  Uh  in  terms  of  Green’s 
functions  to  derive  our  stability  results.  Other  works  for  the  Helmholtz  equation  with 
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radiation  boundary  conditions  of  the  type  (1.2a)  —in  higher  dimensions  too—  include 
[AK],  [G],  [HHl],  [BS]  [IB3],  [BIPS].  Among  these  works,  [AKS],  [HE],  [BS]  and  [BIPS] 
axe  devoted  to  the  analysis  of  non-standard  finite  element  methods  which  have  been 
designed  to  minimize  the  “pollution  effect”  observed  for  the  standard  Galerkin  method 
for  moderate  to  high  wave  numbers.  A  detailed  analysis  of  the  behavior  of  the  Galerkin 
finite  element  approximations  in  the  case  of  high  wave  numbers  is  presented  in  [IBl- 
3]  and  [BS].  For  a  comparison  of  boundary  element  and  finite  element  methods  for 
time-haxmonic  acoustic  wave  equations  cf.  [HH2]. 

2.  Analysis  of  the  Continuous  Problem 

We  first  observe  that  if  is  a  solution  of  (1.6)  then  its  components  satisfy  (1.1-2): 

Lemma  2.1.  Let  J-  £  L^.  Assume  that  (1-6)  has  a  solution  U  =  {p,u,^  G  H.  Then 
U  e  :=  X  H^{^2)  X  and  p,u  and  p  satisfy  (1.1-2),  where  the 

derivatives  in  these  equations  are  derivatives  in  the  sence  of  distributions. 

The  main  part  of  this  section  is  devoted  to  the  proof  of  the  following  stability  result: 

Proposition  2.1.  IfU  is  a  solution  of  problem  (1.6)  for  then 

(2.1a)  ll^^lli.W  <  C*! Iloilo, w, 

where  Ci  is  a  positive  constant  independent  of  k.  Also,  for  .^  £  'W 

(2.1b)  ll^lkw  ^  Ci-\\J^\\i,n , 

We  shall  state  two  consequences  of  this  result  as  Theorems  2.1  and  2.2  below.  The 
proof  of  Proposition  2.1  is  given  after  the  proof  of  Theorem  2.2. 

Theorem  2.1.  Let  T  £  L^.  Then  there  exists  a  unique  solution  of  (1.6).  Furthermore 
14  £  ^ ,  and 

(2.2)  ||Wei||o,-h  <  ^*2(1  +  ^)||'^||o,7i  • 

Proof  The  uniqueness  follows  from  Proposition  2.1.  Now,  using  the  trace  inequality: 

(2.2)  b(a)P  <e||ur||i*(/)  +  ^C'||u|li2(j),  veH\l),e>0 

where  a  is  one  of  the  endpoints  of  the  interval  J,  and  the  arithmetic-geometric  mean 
inequality,  we  easily  verify  that  E(’,  •)  satisfies 

(2.3)  ReB{U,U)>Co\\U\\l.^^--po(l  +  k^)\\U\\l^  WeH, 
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with  positive  constants  Co  ,/xo  independent  of  k.  Therefore  B(-,-)  satisfies  a  Garding- 
type  inequality.  It  is  well  known  that  (2.3)  implies  the  alternative  statement:  Either 
the  problem 

B(i{,v)  =  o  Wen, 

has  a  nontrivial  solution  or  (1.6)  has  a  solution.  Therefore  since  we  have  established 
uniqueness,  (2.3)  implies  existence  of  the  solution  for  the  problem  (1.6).  The  estimate 
(2.2)  follows  from  Lemma  2.1,  Proposition  2.1  and  (1.1).  □ 

Let  II  •  ll'i^^  be  the  dual  norm  of  ||  •  ||i,-h  defined  by  ||V||'i =  sup^g^ 
have  the  following  result: 

Theorem  2.2.  The  bilinear  form  5(-,-)  defined  on  x  satisfies 


(2.4) 

|B(W,V)|<;8||M||i,w||V||,,« 

w,ven, 

and 

(2.5) 

ReS(M,V)^  1„„„ 

sup  >7, 

ve-w  ll^lli.-w  * 

men. 

where  ^  and  7  are  two  positive  constants  independent  of  k.  Furthermore,  if  U  is  the 
solution  of  (1.6),  we  have  the  estimate 

(2.6) 

Proof.  The  inequality  (2.4)  follows  (similarly  as  (2.3)  )  from  (2.2),  the  arithmetic- 
geometric  mean  and  Cauchy-Schwarz  inequalities.  For  the  proof  of  (2.5),  let  14  be  a 
given  element  of  H.UV  =  where  Z  eU  shall  be  appropriately  chosen  bellow, 

we  have 

B{U,V)  =  ^B{U,l4)+BiU,Z) 

1  /  1  +  P  \ 

=  ^  {B{U,U)  +  fioil  +  k^mWln)  +  [B{U,Z)  - 

=:I+II. 

We  choose  Z  such  that  the  part  II  vanishes.  I.e.,  Z  shall  be  the  solution  of 
(2.7)  B(^,  Z)  =  no^  ^^^—{^,U)o,n 

(2.7)  is  a  dual  problem  to  (1.6)  and  one  can  verify  that  the  results  of  Proposition  2.1 
and  Theorem  2.1  hold  for  this  problem  as  well.  Therefore  Z  is  uniquely  determined  as 
solution  of  (2.7),  and  furthermore  satisfies  the  estimate  (since  U  eH) 

(2.7)  ||Z||,,«<C.i,i.i^||M||i,«. 
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With  this  choice  of  Z  we  have 

(2.8)  ReB{U,V)>^C4U\\ln. 

where  V  G  H  satisfies 

||V||i,-w  <  -j^\\U\\i,n  +  ll-^lli.w 


Therefore  (2.8)  implies 


ReB(W,V)  >-^^^\\U\\ln 

>Co||ZY||i,-w^^||V||i,-w  =  ’ 

and  (2.5)  follows.  The  estimate  (2.6)  is  a  consequence  of  [BA,  Theorem  5.2.1]  and  of 
(2.4),  (2.5).  n 

We  now  prove  Proposition  2.1: 

Proof  of  Proposition  2.1.  We  write  B{U,  V)  =  5i(ZY,  V)  +  £2(^5  V)  +  .B2(W,  V),  where 

Bx{U,  V)  =  /  PtIx^x  I  Pldx  -  Pu{xi)q(xi)  -  ikp{0)q{0) , 

Joi  Jill 


(2.9)  B2{U,V)=k^  f 

Jii 


aUxV 


^dx  —  k*  guvdx  —  k“^ p{xi)v{xi)  +  k^p{x2)v{x2)  , 

Jil2 

-  k^  I  pqdx  +  k^u{x2)q{x2)  -  ikp{L)q{L) . 


BziU,  V)=  f  Pxqxdx  -k^  I  pqdx  +  k^u{x2)q{x2)  -  ikp{L)q{L) . 

In  the  sequel  we  assume  that  if  is  a  solution  of  (1.6)  and  T  G  I?.  By  Lemma  2.1  we 
have  U  and  the  components  of  U  satisfy  (1.1-2).  Therefore  the  function  V  with 

components 

q-{x-  xi)px  , 

(2.10)  u  =  (a:  -  X2)ux  , 

q  =  {x-  X2)px  , 

is  an  element  of  H.  We  shall  evaluate  ReB(if,  V).  We  observe  first  that  (?(xi)  =  0  ): 

Bi(if,  V)  =  /  Pxqxdx  I  Pqdx  -  Pu{xi)q{xi)  -  ikp{0)q{0)  , 

Jai  Jill 

=  f  px{x-xi)pxxdx+  I  \Px\‘^-k^f  p{x-xi)pxdx  +  ikp{0)xipx{0); 
Jqx  JQi 
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hence,  using  (1.2), 


ReBi{U,V)=  f  \px\^dx  +  Rel  f  {x  -  xi){\p^\%dx 
-Re^P  f  ix-xr){\p\^)rdx-xMO)\^ 

2  Jill 

=  J  +  \Pxfdx  +  xi\pj;{0)\^^ 
-Re^fc^j-  J  bprfa:  +  a:ilp(0)p| -a:i|px(0)|^ 

=-  /  \Px\'^dx  +  ^P  f  bprfa:  -  ixi|px(0)p  -  ^fc^xi|p(0)P  . 
2  Jdx  ^  JQt  ^  ^ 


Since  |pi(0)l  =  fcb(0)|  we  conclude 


(2.11)  ReBi(U,V)  =  l  f  \Px\^dx  +  h^  f  [pl'^dx  -  xik^\piO)\^  . 

2  Jill  2  Jill 


Similarly, 


B2{U,V)  /  auxVxdx  -k‘^  I  Quvdx  -  k^p{xi)v{xi)  +  Pp{x2)v{x2) , 

Jn^  J^2 

— /  dUxi^x  X2^'dxxdx  “f"  k  j  ci|ux I  dx 
Jui  dil2 

—  k^  I  gu{x  —  X2)uxdx  —  k^p{xi){xi  —  X2)ux{xi)] 

JQ2 


hence,  using  (1.2), 


ReB2iU,V)=k^  f  a\ux\'^dx  +  Re]-k^  f  a{x  -  X2){\ux\^)xdx 

Jq2  ^  dU2 

-Re-A:^  /  g{x  —  X2){\u^)xdx -Rek'^p{xi){xi  —  X2)ux{xi) 

2  Ju2 

=\k"^  f  a\ux\‘^dx +  ^Pa{x2  -  xi)\ux{xi)\^ 

2  Ja2  ^ 

+  f  g\u\^dx  -  ^k*g{x2  -  a:i)|u(a:i)p  +  Pa{xi  -  a;2)lux(3;i)h 

2  in,  2 


Therefore 


ReR2(^,V)  =^A:^  f  aluxl^dx  +  ^A:^  /  g\ufdx 

2  Jn2  ^  J^2 


^k'^gix2  -  xi)\u{xi)\^  -  ^Pa{x2  -xi)|ux(xi)p  . 


(2.12) 


It  remains  to  evaluate  BziJA^V)  : 

Bz{l{,V)=  f  Mx-^2)Pxx<i^+ I  /  p(a;-a;2)Pxt^a;-iA:p(L)(L-X2)Px('^) 

Jqi  JQs  Jilz 

and,  using  (1.2), 

Re53(ZY,V)=y  \Px\^dx  +  Re^l^-  \Px\^dx  +  (L  -  X2)\ML)\^^ 

-  Re^A:^  |-  J  {^^dx  +  (L  -  X2)|p(I')p|  -  (L  -  X2)\px{L)\^ 

=  \  j  \Px?dx  +  ]-k'^  j 

2  Jn,  ^  J^3 


fQ3 

-l(L-  »2)|p.(I)P  -  lk^(L  -  . 


or,  since  \px{L)\  =  A:|p(X)i, 

(2.13)  RcBi{U,V)  =  \  j  \Px\^dx^-\k‘^  f  dx  -  {L  -  X2)k'^  \p{L)\^  ■ 

2  y«3  2  7^3 

From  equations  (2.11)  and  (2.13)  it  is  clear  that  we  need  to  estimate  |p(0)|,  \^L)\. 
To  this  end,  we  observe  that 

(2.14)  lrxiB{U,U)  =  k\p(0)\^  +  k\p{L)\^  , 
and  therefore 

(2.15)  i"|p(0)|=  +  *:"|p(i)l"  <  k\(y^,U)n.n\ ■ 


Now,  let  V  =  Vi  +  V2  +  Vs  where 

Vi  =  (9,0,0),  V2  =  (0,u,0),  V3  =  (0,0,^ 

{q,v  and  9  are  given  by  (2.10)).  The  relations  (2.11)  and  (2.13)  are  still  valid  if  we 
replace  V  by  Vi  and  V3,  respectively.  Hence,  in  view  of  (2.15),  we  obtain 


(2.16) 


^\\pxf  +  ^^^\\pf  ^  |(^,Vi)o,-h|  +a:iA:|(/’,W)o,w|, 

^llP^ir  +  ^3)0, Til  +{L-  X2)k\{T \  • 


We  turn  now  to  (2.12).  In  order  to  estimate  the  terms  at  aii,  we  rewrite  them  using 
the  transmission  condition  (1.2b)  as 


(2.17) 


\k'^g{x2  -  xi)|u(xi)p+-Pa(a:2  -  Xi)|«x(xi)l^ 


=  ^g{x2  -  Xl)lpx(xi)r  +  -  Xl)|p(xi)P  . 
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To  bound  the  right-hand  side  of  this  equation,  we  proceed  as  follows:  Let  q'  =  xpx  6 
and  =  (9',0,0).  Using  the  relation  k‘^u{xi)  =  Px{x\)  and  following  the  steps 
of  the  proof  of  (2.11)  we  conclude 


ReBi{U,V[)=  \Px\^dx+Re^  M  \Px\  dx  Xl|pi;(xi)|  ^ 

-ReiA;^|-  J  \p\^dx  +  \p{xi)\‘^^  -RePxiu{xi)p^{xi) 

=1  f  \Pxfdx  +  \p  f  \p\^dx-^xi\pxixi)\^ -^Pxi\p{xi)f  . 

2  Joi  2  2  ^ 


Hence,  using  (2.15),  we  obtain 

Wxi)|“  +  t=|p(xi)|=  <i  l|p.r +  i  + 


(2.18) 


XI 


XI 


<1^  +  |(^,  VOo,7i|)  +  2k\iT,U)o,n\ ■ 

ZiX 


Summarizing,  (2.12),  (2.17)  and  (2.18)  give  the  estimate 


+  if‘ll«||"  <|(:^,V2)o,«|  +  C.|(;f,V,)„,«| 

+  C,\(r,V[)o,-H\  +  C.fc|(:F,M)o,«| , 


where  C*  is  a  constant  (not  necessarily  the  same  in  two  different  places)  independent 
of  k.  Recalling  the  definitions  of  Vj,  VJ  we  have 

(2.19)  k^UxW^  +  k^uf  <  a{P\\f\\  llu.ll  -f  Ikill  Ib.ll)  +  ak\{:F,K)o,-H\  > 

where 


(2.20) 


k\(T,U),.n\  <*:|ki||  IIpII  +  <:’ll/ll  ll>‘ll  +  k^g4  liai 

<c.(||s,f  +  /fc^ll/ll  +  IIpjII)  +  £(t"l|p||  +  fc‘ll«ll  +  k-^m) 


Handling  analogously  the  terms  of  (2.16)  and  sumnaing  we  conclude: 


l|p.r  +  *"l|pf  +  k'‘\\n,f  +  +  llp.lp  + 

<C.(|kf +  t"ll/ll  +  IIS2ll), 


or 

ll^lll.W  —  11*^110, 5 

which  is  the  relation  (2.1a).  If  ^  we  have  l|.F||o,?<  <  il|.^||i,7i  and  (2.1b)  foUows.  □ 

Remark.  It  is  clear  from  the  proof  of  Proposition  2.1  (cf.  (2.12),  (2.17-9)  )  that 
the  constant  Ci  of  (2.1a, b)  tends  to  infinity  as  a  — »•  0,  in  particular  Ci  =  O(^). 
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2.  Finite  Element  Approximations 

In  this  section  we  shall  consider  numerical  approximations  to  the  solution  of  problem 
(1.6).  If  Si,  i  =  1,2,3,  are  three  finite  dimensional  subspaces  of  i  =  1,2,3, 

and  Sh  is  the  space  5ft  =  5^  x  5^  x  5|  we  have  5ft  C  H.  We  shall  seek  approximation 
£  5ft  to  ZY  defined  as  the  solution  of  the  problem:  Find  Uh  €  5ft  such  that 


(3.1)  B{Uh,^)  =  V<fG5ft, 

where  B{-,  •)  is  defined  in  (1.6b). 

The  finite  element  spaces.  Our  basic  assumption  on  the  finite  element  spaces  5^, 
i  =  1,2,3,  is  that  they  consist  of  continuous  piecewise  polynomials  of  degree  s,  s  > 
So  >  1.  These  functions  are  constructed  using  a  fixed  but  arbritrary  partition  of  Ui  into 
intervals  {Pj}j.  Note  that  the  degree  of  the  polynomials  may  vary  in  the  intervals 
Let  h’  =  maxj  |lj|  and  h  =  max{/l^/l^/l®}.  A  consequence  of  standart  approximation 
properties  for  these  spaces  is  the  relation:  For  any  V  G  there  is  an  element  ^  G  5ft 
such  that 

u2-m 

(3.2)  <  C - ||V..||o,w  ,  m  =  0, 1 , 

^  'So 


where  =  V  and  =  Vx  •  Of  course  (3.2)  gives  a  crude  bound  for  inf,i^g5^  ||V  - 
but  we  use  this  estimate  as  a  tool  in  the  proofs  below.  In  particular,  we  first 
show  that  the  bilinear  form  B{‘,  •)  satisfies  a  discrete  analog  of  the  BB  condition  (1.10) 
(Theorem  3.1),  and  also  that  two  quasioptimal  estimates  of  the  type  (1.12)  axe  satisfied 
((3.9)  and  Theorem  3.2).  An  essential  tool  in  the  proofs  of  these  results  is  the  stability 
estimate  of  Proposition  2.1. 

Theorem  3.1.  Assume  that  h  is  small  enough  such  that  <  c,  where  c  is  an 

appropriately  chosen  constant.  Then  there  exists  a  constant  7  >  0,  independent  of  h 
and  k,  such  that 


(3.3) 


Re  B{Uh,Vh)  ^  II 

— Nvli - >  7Tll^ft||l,'N, 


VZ/ft  G  5ft . 


Proof  We  first  consider  the  auxiliary  bilinear  form 

(3.4)  A(ZY,V)  =  B(ZZ,V)+//o(l  +  A:^)(^,V)o,w,  U,V£H, 

where  /xq  is  the  constant  of  (2.3).  It  is  easy  to  see  that  there  is  a  constant  (3 a  >  0, 
independent  of  such  that 


\A{U,V)\<^A\\Uh,n\\V\\i,n. 
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Also  (2.3)  implies  that  A{-,  •)  is  coersive  with  respect  to  ||  •  Therefore  the  projection 
Px  •  W  is  well  defined  by 


(3.5)  A{^,  PaV)  =  A{^,  V)  WeSh. 

For  the  proof  of  (3.3),  let  Uh  be  a  given  element  of  Sh-  Let  Vh  =  +  PaZ,  where 

Z  e  shall  be  defined  below.  For  this  choice  of  Vh  we  have 

(3.6) 

B{Uh,Vh)  =  -^BiUhM)  +  B{U,PaZ) 

=  p  [BiUhiUh)  +  ^o(l  +  ^^)||^/i||o,'h)  +  ^B{Uh,PAZ)  —  fio  p  II^aIIo.w^ 

=  p  +  /*«(!  +  *^*)II«).|I2,«)  +  {b(U^,  Z)  - 

+  /io(l  +  Z  —  PaZ)(s,'h  • 

Choosing  now  Z  ^  'ht  io  be  the  solution  of 

(3.6’)  =  V<f€W, 

the  second  term  of  the  last  equality  of  (3.6)  vanishes.  (As  we  have  seen  in  the  proof  of 
Theorem  2.1  Z  is  uniquelly  determined  as  solution  of  (3.6*),  and  furthermore  satisfies 
the  estimate  \\Z\\i,h  <  Ci\^iQ^^\\Uh\\i,n  •)  Therefore  it  remains  to  estimate  the  term 

{llh,Z  -  PaZ)o,h- 

We  shall  use  the  following  lemma,  the  proof  of  which  is  given  below. 

Lemma.  3.1.  If  Z  ^  Ti  and  PaZ  is  defined  by  (3.5),  then  we  have  the  estimate 

\\Z-PAZ\\o,n<C^\\Z\U,-H, 

where  the  constant  C  >  0  is  independent  of  h  and  k. 

Now,  since  \\Z\\r,-H  <  C\\\Uh\\i,H  and  \\Uhh,-H  <  ,  we  obtain 

/xo(l  +  k^)\{Uh,Z  -  PAZ\n\  <C{\  +  k^)\\Uh\\o,n^\\Z\\x,-H 

So  ’  So 


Therefore,  if  k^-^  is  small  enough,  we  obtain 


(3.7) 


IUB(%,Vk)  >pC„||W||f,« 


C^IIMII 

So 


2 

i.n 


12 


where  V/i  G  <5/1  satisfies 


||V|l||i,J(  <  p|iw»||i,7<  +  WPaZWi.-h  <  p||%||i,H  +  C||Z||i,7i 

<^(Ck  +  i)  ||W»||,,« 

<^Czk\\U4^,n. 

Therefore,  (3.7)  implies 

which  yields  (3.3).  □ 

Proof  of  Lemma  3.1.  We  shall  use  a  classical  duality  argument.  Let  p  :=  Z  —  PaZ  and 
3^  G  Tf  be  the  solution  of 

(3.8)  A{^,  #)  =  k^ip,  ^)o,«,  Wen. 

This  problem  has  a  unique  solution  with  the  regularity  3^  G  n^ .  Also,  using  the  coersivity 
of  A{-,  •)  and  equation  (1.1),  we  easily  conclude  that 

Ill'll Ho.-H  <  ^^||p||o,-H- 

Therefore,  using  (3.2)  for  a  properly  chosen  iP*  €  Sk. 

<  ~  ^h\\i,n\\Ph,'>i  ^  C—\\^xx\\o,n\\ph,'H 

So 

<CP^\\p\\oMi,n; 

So 


hence, 

WI.,«<cflWI..«<Cp||z||..-w, 

5o  5o 

since  from  the  definition  of  Pa  we  have  ||Pa-2^||i,'W  <  .  □ 

An  application  of  Theorems  2.1  and  3.1  and  of  [BA,  Theorem  6.2.1]  gives  the  follow¬ 
ing:  If  the  hypotheses  of  Theorem  3.1  hold,  then  Uh  is  uniquely  determined  as  solution 
of  problem  (3.1)  and 

(3.9)  WU  -  llh\\i,-H  <  Ck  inf  \\U  -  <?||i,-w  ,  C>0  independent  of  k. 
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where  tl  is  the  solution  of  (1-6). 

The  estimate  (3.9)  is  a  quasioptimal  result  of  the  finite  element  approximation  with 
respect  to  ||  •  ||  i  .  It  is  possible  however  to  improve  this  resxilt  in  terms  of  the  dependence 
of  the  boimd  on  k.  Indeed,  in  the  proof  of  the  following  result  we  shall  replace  the 
constant  Ck  in  (3.9)  by  a  constant  independent  of  k.  A  comparison  of  (3.9)  and  (3.10) 
below  indicates  that  (3.9)  might  be  valid  under  a  weaker  restriction  than  Ph  <  c,  c 
small.  This  assumption  is  consistent  with  avaihable  results  for  the  single  reduced  wave 
equation,  cf.  [IB1,2],  where  the  analog  of  (3.3)  and  therefore  of  (3.9)  is  proved  (in  a 
uniform  mesh)  under  the  weaker  condition  that  kh  remains  bounded.  The  proof  of  the 
Theorem  3.2  is  based  on  a  argument  due  to  Schatz,  [S].  The  use  of  the  stability  estimates 
of  Section  2  is  essential  here. 

Theorem  3.2.  Assume  that  the  hypotheses  of  Theorem  3.1  hold.  Then 

(3.10)  WU  -  lfh\\i,n  <  C*  inf  \\U  -  ,  C*  >  0  independent  of  k, 

where  lA  and  lAh  the  solutions  of  (1*6)  und  (3.1),  respectively. 

Proof.  Let  ^  ^7{  he  the  solution  of  the  problem  {e=U  —  Uh) 

(3.11) 

As  we  have  seen  in  Section  2,  this  problem  has  a  unique  solution.  Now,  for  any  Th  €  Sh^ 

The  results  of  Theorem  2.1  hold  for  'P  as  solution  of  an  adjoint  problem  to  (1.6).  There¬ 
fore  P  G  TO"  and 

\\^rx\\o,-H<C{l  +  k)P\\e\\o,n- 

Hence,  by  (3.2), 

•So 

<Ct'-(l  +  t)||e||o,«|le||,,„. 

^0 

or 

(3.12)  ||e|lo,«  <  C'(l-hfc)^||el|i,w. 

Since  the  bilinear  form  B(-,  •)  satisfies  (2.3),  we  have,  for  any  P  e  Sh  , 

Co\\e\\ln  =  B{e,e)  +  fxo{l  +  k^mU 

=  B{e,U  —  <f)  +  fio{l  +  A:^)|le||o 

P 

<  /^lkl|i,?i||W  -  +  /^o(l  +  ^^)(1  +  ) 
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4.  Numerical  Evaluation 


To  demonstrate  the  principal  numerical  effects  of  Galerkin  FE-solutions  to  the  Helmholtz 
equation,  we  first  consider  a  related  uncoupled  problem,  [IBl], 

u"  +  k^u  =  1  on  fii  =  (0, 1) 

(4.1)  «(0)  =  0 

u'(l)  —  iku{l)  =  0. 

For  the  finite  element  solution,  fl  is  covered  by  uniform  mesh  Xh  of  n  elements  with 
stepsize  h  —  n~^.  The  subspace  of  piecewise  polynomial  functions  of  uniform  degree 
p  with  nodes  in  Xh  is  called  5'^.  The  error  is  measured  in  if^-seminorm  |  •  |i  which, 
for  this  problem,  is  equivalent  to  the  full  /f^-norm.  The  optimal  convergence  rate  is 
then  displayed  by  the  -projection  (in  the  seminorm)  of  the  exact  solution  onto  5^. 

It  is  well  known  that  this  best  approximation  interpolates  the  exact  solution  on  Xh- 
In  Fig.  1,  the  errors  of  best  approximation  and  finite  element  solution  are  shown  for 
k  =  2,k  =  10  and  k  =  100. 

We  observe  that  for  all  k  the  Galerkin  FE-solution  reaches  optimal  convergence  if  h 
becomes  small.  However,  for  medium  and  large  k,  the  range  of  optimal  convergence  is 
preceded  by  a  range  in  which  the  FEJ-solution  is  numerically  polluted,  cf.  [BS],  [IB3].  In 
this  range,  the  finite  element  error  differs  significantly  from  the  optimal  approximation. 

In  particular,  it  was  shown  by  numerical  experiment  that  a  constraint  of  the  form 
h?k  =  const  is  indeed  necessary  for  quasioptimal  behaviour  of  the  FE-solution  if  So  =  1? 
cf.  [IBl]  -  see  Fig.  2. 

For  sufficiently  regular  solution,  oscillating  with  frequency  fc,  the  relative  error  in 
If ^ -norm  can  be  estimated,  generally,  as  [IB2] 

ei  <  CiihkY”  +  C2k{khf^\ 

It  was  shown  in  numerical  experiments  that  this  estimate  is  sharp  in  the  preasymptotic 
range  of  convergence.  Hence,  raising  the  order  of  polynomial  approximation  So  with 
hk/so  =  const,  the  gain  in  performance  of  the  Galerkin  FEM  for  numerical  solution 
of  the  Helmholtz  equation  is  twofold:  increase  of  the  order  of  asymptotic  convergence 
and  decrease  of  the  pollution  in  the  preasymptotic  range.  This  is  illustrated  by  the 
computational  results  shown  in  Fig.  3.  The  vertical  line  at  the  relative  error  ei  =  0.1 
indicates  the  decrease  in  the  number  of  elements,  needed  for  a  given  accuracy  as  the 
degree  of  approximation  is  increased.  For  a  more  detailed  discussion  of  the  numerical 
effort  involved,  see  [IB2].  The  “bumps”  in  the  plots  are  an  effect  of  local  approximation 
if  the  meshsize  coincides  with  the  wavelength,  [IB2]. 

Turn  now  to  computational  residts  for  the  model  problem  (1.1-2)  of  fluid-structure 
interaction.  The  numerical  evaluation  has  been  carried  out  with  the  data: 

•  The  domains  are  fli  =  (0,3);  CI2  =  (3,6);  fls  =  (6,9). 
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•  In  the  first  region,  a  source  is  present  in  the  form  of  a  step-function 

f  0  in  (0, 1) 

^1  =  <  1  in  [1,2] 

[o  in  (2,3) 

•  No  load  is  given  in  the  sohd  and  in  the  ‘right’  fluid  region: 

f  =  92  =  0. 

All  domains  are  covered  by  uniform  (in  h  and  s)  meshj  the  computations  are  carried 
out  with  10, 20, 40, ... ,  320  elements  per  region.  The  degree  d  of  polynomial  approxi¬ 
mation  was  varied  from  1  to  7.  The  errors  are  measured 

1.  pointwise  as  absolute  values  in  xi  =  1.5,  X2  =  4.5  and  xz  =  7.5. 

2.  In  the  ifi-norm  on  the  first  region. 

In  the  second  case,  the  FE-error  was  compared  to  the  error  of  the  iT^-projection.  A 
review  of  the  plots  for  all  computations  shows  that  the  coupled  problem  displays  the 
same  principal  effects  as  the  uncoupled  problem: 

1.  The  error  is  smaller  than  100  %  only  if  the  meshsize  exceeds  a  critical  number 
which  grows  over-proportionally  with  the  wavenumber  k  —  see  Fig.  4. 

2.  On  coarse  mesh  and  for  medium  or  large  k,  the  pollution  term  dominates  the  FE- 
error.  For  small  h,  the  FE-solution  is  quasioptimal  with  an  optimality  constant 
not  depending  on  k  (Theorem  3.2).  This  behaviour  is  shown  in  Fig.  5  where 
the  relative  errors  of  the  -projection  and  the  FE-solution  are  compared  for 
different  k  and  different  degrees  of  approximation.  The  quasioptimal  behaviour  of 
the  solution  for  sufficiently  small  ^  can  be  seen  in  Fig.  6. 

3.  Increasing  the  polynomial  degree  So,  one  decreases  the  influence  of  the  pollution 
term  on  the  error  and  increases  the  rate  of  convergence  -  see  Fig.  7. 

The  graphics  shown  here  for  demonstration  are  representative  for  the  whole  series 
of  computationeil  results.  In  particular,  no  significant  dependence  of  the  error  on  the 
region  Hi  was  observed. 
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Relative  error 


Figure  1:  Errors  of  /f^-projection  versus  error  of  finite  element  solution  for  Dirichlet 
problem;  error  in  ff^-seminorm;  wavenumbers  =  10,  A:  =  50  and  k  =  100. 
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Figure  2:  Constraints  of  the  form  k^h  =  a  for  a  =  1  and  a  =  0.1  in  error  plots  for 
k  =  10,50,100  and  k  =  200. 
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Relative  error 


Figure  3:  Relative  error  of  the  FE-solution  versus  error  of  ff^-projection  for 
A:  =  50  and  So  =  1, 2, ...  6. 
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number  of  meshpoints 


Figure  4:  Relative  error  of  the  finite  element  solution  e 
Tiex(®o)|/ke*(®o)|>  for  the  coupled  problem  at  Xj,  =  1 

10,40,100,400.  (a)5„  =  l. 


=  \UfeiXo 

.5,  for  k 
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Figure  5:  Pollution  effect;  Relative  error  of  ff*-projection  and  finite  element 
solution  for  the  coupled  problem;  error  in  if^-norm  for  k  =  10, 30, 100  and 
jfc  =  500.  (a)  jb  =  10,5<,  =  2 


24 


■  1  t  ■ 


10 

Number  of  elements 


10 


30, io  =  2 


25 


Figure  5;  (c)  k  =  100, 5o  =  4 
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Figure  5:  (d)  k  =  500,  So  =  5 


27 


Number  of  elements 


Figure  7:  p-  convergence:  Relative  error  in  H^-norm:  error  of  the  finite 
element  solution  vs.  error  of  best  approximation  for  k  =  300  and  So  =  1 . . .  7 
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